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1.0 


BACKGROUND AND INTRODUCTION 


Since the inception of this project, many existing mathematical 
models of major physiological subsystems have been identified and 
simulation models have been developed, modified, and validated. 

Several new models have also been developed and tested in this pro- 
gram. This contract encompasses the first version of the whole- 
body algorithm which integrates these subsystem models to study major 
physiological Interactions. Recent modeling activity has been to 
assemble, develop, and Integrate physiological models and to collect 
pertinent literature from the last few years. One natural evolution 
from this point is the consideration of possible applications of simula- 
tion models in diagnosis and therafieutics. 

This report presents the highlights of current development of 
models for clinical applications and results of related research studies. 
Since the major emphasis of this contract has been development, assem- 
bly, and integration of simulation models of various physiological sub- 
systems, the application of these models for diagnostic and therapeutic 
problems is in a very early stage. A few clinical areas, however, have 
been studied and are discussed in the following sections. 
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2.0 CURRENT APPUCATIONS OF SIMULATION MODELS FOR 
CLINICAL RESEARCH 

In reviewing applicable literature, it is clear that various physio- 
pathological rtates have been studied by many Investigators, but their 
selection of subjects, healthy or 111, are random and, of course, are 
limited to the specific regulation/control responses in one or a few 
body subsystems of Interest to the invectigatlon. The application 
of models to the circulatory system was selected here. 

2.1 CIRCULATORY MODELS FOR CUNICAL RESEARCH 

2.1.1 Tilt Model Simulation of Orthostatic Intolerance with Hemorrhage 

Despite extensive research, the causes of many common distur- 
bances of the human circulation remain obscure. The events leading 
to changes in blood pressure are not completely understood even for 
simple changes In posture. Recent space flights of long duration 
causing extended exposure to a weightless environment, and the 
resultant orthostatic fall in blood pressure upon reentry have renewed 
interest in gaining a better understanding of the underlying mechanisms 
of this phenomenon. The same orthostatic Intolerance has been 
observed on first standing after prolonged bed rest. Computer simu- 
lation of the circulation and its controls permits a more rigorous 
analysis of these mechanisms. 

The force of gravity causes blood to pool in the lower body (Gauer 
(1963(b)) which causes a transient fall in pressure to be sensed at the 
baroreceptors in the carotid artery and the aortic arch. When the 
pressure falls at the baroreceptor, the sensory nerves will deci*ease 
firing and there will be Increased stimulation of the sympathetic fibers 
with decreased stimulation of the vagal fibers. This autonomic activity 
will produce vasoconstriction. Increased heart rate, an Increase in 
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stroke volume, and an increase in peripheral resistance (see Figure 
1). 'Hie increased arteriolar resistance acts to restore arterial 
pressure as discussed by Henry (1955) and thus, a new steady-state 
is reached. SJostrand (1952) Indicates that stroke volume is depen- 
dent on cardiac filling, which responds to the decrease in central blood 
volume because of the legward :«hlft due to gravity, and to a lesser 
extent on autonomic activity. ‘Hie stroke volume is decreased due to 
the reduced central blood volume more than the heart rate is increased 
by the autonomies and the result is a decrease in blood pressure. It 
can be seen that without the autonomic reflex actions to restore cardiac 
output and increase peripheral resistance and thus restore arterial 
pressure through the baroreceptor feedback loop, the system would 
become unstable and one would probably faint every time the upright 
posture was assumed. All of these essential compensations depend 
on actions of the autonomic nervous system, which responds to the 
depletion of central blood volume and to a transient fall in arterial 
pressure. The model's response to a change in posture from supine 
to upright (see Figures 2 through 4 for normal blood volume) resembles 
the behavior of the human system as shown in Table 1. 

It has been observed that there is a moderate decrease in total 
blood volume after extended exposure to a weightless environment or 
prolonged bed rest (Berry (1967), Gauer (1963 ), and Waterfield (1931) . 
This may be due to Increased capillary filtration out of the central blood 
volume due to the increased pre-capillary pressure as a result of the 
headward migration of blood which occurs with the removal of the 
gravity force in the lower body and an observed reduction in red cell 
mass. The orthostatic Intolerance which results after such expos- 
ure can then be simulated by a reduction in total blood volume (slow 
hemorrhage). The model's response to an upright tilt after a one 




FIGURE 1. BLOCK DIAGRAM REPRE'JENTING BARORECEPTOR 

AND AUTONOMIC RESPONSE TO A CHANGE IN 
BLOOD PRESSURE 











TABLE 1. TILT MODEL SIMULATION RESULTS COMPARED T'' 
EXPERIMENTAL DATA FOR PASSIVE UPRIGHT 
TILT BEFORE AND AFTER ONE LITER HEMORRHAGE 



SUPINE-NORMAL B.V. 

CHANGE AFTER 90° 
TILT (NOR. B.V.) 

SUPINE - AFTER 
ONE LITER HEMORRHAGE 


MODEL 

RESULTS 

EXPERIMENT 

DATA* 

MODEL 

RESULTS 

EXPERIMENT 

DATA* 

MODEL 

RESULTS 

EXPERIMENT 

DATA* 

MEAN ARTERIAL PRESSURE (mm Hg) 

87 

100 

+5.0 

-5 TO +5 

79 

85 TO 100 
• 

CARDIAC OUTPUT (liters/min) 

6.6 

5.86 

- .4 

- .86 

5.5 

4.5 TO 6.0 

HEART RATE (bcats/min) 

63 

64 

+ 12.5 

+16.0 

71 

65 TO 90 

STROKE VOLUME (ml) 

106 

92 

-23 

-28 

78 

60 TO 90 

LOWER BODY BLOOD VOL (ml) 

500** 

1400 

+380** 

+700 

470 

- 

TOTAL BLOOD VOLUME (liters) 

5.0 

5.7 

0 

0 

4.0 

4.7 


* EXPERIMENTAL DATA FROM REFERENCES Gauer (1963) Weissler (T957), Wriaht (1965), 

Shenkin (1944), and W-rren (194 5T. 

MODEL LOWER BODY VOLUME INCLUDES ONLY THE LEGS 


HtMtT **TE WATS VM) LCC tt.000 VOC (UTEISl 



FIGURE 2. TILT MODEL OUTPUT OF LEG BLOOD 

VOLUME AND HEART RATE FOR 
UPRIGHT TILT SIMULATION BEFORE 
AND AFTER ONE LITER HEMORRHAGE 




FIGURE 3. TILT MODEL OUTPUT OF CARDIAC 

OUTPUT AND STROKE VOLUME FOR 
UPRIGHT TILT SIMULATION BEFORE 
AND AFTER ONE LITER HEMORRHAGE 





FIGURE 4 . TILT MODEL MODEL OUTPUT OF BLOOD 

PRESSURE AND MEAN ARTERIAL 
PRESSURE FOR UPRIGHT TILT 
SIMULATION BEFORE AND AFTER ONE 
LITER HEMORRHAGE 
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liter blood loeo If compared to the normal blood volume caae In 
Flgurea 2 through 4. After removal of one liter of blood, the central 
and peripheral compartmenta are depleted proportional to their com- 
pllancea. In the aupine poaltlon (compare pre-tllt values) after this 
blood loss, the autonomic relfexea are able to maintain arterial pres- 
sure while venous return and cardiac output are moderately reduced 
(see Table 1). The system shows good stability so long as the recum- 
bent posture Is maintained. U the subject attempts to stand after the 
loss of one liter of blood, central blood volume becomes critically 
depleted (note the drastic reduction in stroke volume) and the blood 
pres^i.re becomes dangerously low. This condition is often more 
than can be compensated for by the autonomies, and vasodepressor 
syncope (V'eissler (1957)) results. These results agree with the 
observations of Shenkln (1944) In men who lost one liter of blood. 

Half of the subjects maintained a stable cardiac output and arterial 
pressure after loss of one liter of blood. Other subjects showed 
effects ranging from quickening pulse to premonitory symptoms of 
fainting. 

The model can also be used to examine more complex circula- 
tory reactions to stress and various changes in autonomic control 
(such as changes in contractility of the heart, increased baroreceptor 
sensitivity, increased peripheral resistance response, changes in 
venous tone, and combinations of the above). 

2.1.2 Long Term Circulatory Regulation 

Most of the basic circulatory models used in experimental appli- 
cations have utilized a basic closed-loop flow system with no leaks in 
order to describe various facets of circulatory function. This kind 
of approach can be Justified when the challenge being simulated is acute 
(1-15 min.), but when longer duration simulations are required, a large 
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number of other regulatory mechanisms must be included in order to 
describe overall circulatory control, even in a crude way. An initial 
attempt to model the major components involved in long term circula- 
tory regulation has been made by Guyton, et al(1972), and this analysis 
has proved to be so useful in such a variety of different situations that 
the significant features of this long term model will be presented here. 

Detailed flow charts and model explanations are available (Guyton 
et al (1972), Guyton et al (1973), White (1973). WTiite (1974)) and will 
not be repeated. The systemic circulation is represented by the 
five-compartment model shown in Figure 5. In this figure, the 
abbreviations used are arterial volume (VAS); venous volume (VVS); 
right atrial volume (VRA); pulmonary arterial volume (VPA); pulmon- 
ary venous and left atrial volume (VLA); pulmonary vascular resis- 
tance (RPT); non-muscle, non-renal vascular resistance (RSN); 
muscle vascular resistance (RSM); renal resistance (RR); and large 
vein resistance (RVG). This model is similar to that pictured in 
Figure 6 except that parallel capacitive beds are not considered, but 
parallel flows are included. The model utilizes similar basic cardiac 
function curves as well, but includes the effects of autonomic stimulation, 
arterial pressure afterload, and cardiac hypertrophy or degeneration on 
the pumping ability of the heart. The unstressed volumes of each capac- 
itive region are controlled by the level of angiotensin D in the blood, the 
level of autonomic stimulation, and by the pressure in the region (through 
stress relaxation). The flow- resistances are controlled by a combination 
of local effects and hormonal effects. The oxygen transport features of 
the circulation are present and hematocrit and red cell control are con- 
sidered. This circulation is not closed, but "leaks" through the capil- 
laries, "excretes" through the kidneys, and "drinks" directly into the 
biood. The blood, composed of plasma (with dissolved proteins and 











FIGURE 6. FLUID C( 

MODEI 
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electrolytes) and red cells, thus serves as a filterable fluid and Figure C 
illustrates the fluid containing reservoirs of the model, llie symbols on 
this figure represent plasma volume (VP), red cell volume (VRC), 
interstitial gel volume (VG), free interstitial fluid volume (VIF), cell 
fluid volume (VIC), pulmonary fluid volume (VPF), water Intake rate 
(TVD), urinary output rate (VUD), red cell production rate (RCl), red 
cell destruction rate (RC2), rate of capillary filtratiou (VTC), lymph 
flow rate (VTL), and the rate of fluid flow into the cells (VID). The 
capillary filtration rate is determined from a whole-body version of 
Starling's relationship which states that net filtration pressure Is equal 
to capillary pressure, plus tissue colloid osmotic pressure, minus inter- 
stitial fluid pressure, and minus plasma colloid osmotic pressure. Pro- 
tein is produced and lost by the body and is distributed between the inter- 
stitial space and the plasma. 

Although the above description is far from complete, even for the 
present rather crude model, it illustrates the variety of component systems 
necessary for long term circulatory control. In fact, at least eight major 
pressure control mechanisms are present in this model. These are the 
baroreceptor mechanism, the chemoreceptor mechanism, the central 
nervous system Ischemic response mechanism, the mechanism of stress- 
relaxation, the capillary fluid shift mechanism, the renin-angiotensin 
mechanism, the renal-body fluid control mechanism, and the aldosterone 
control mechanism. One way of characterizing the effects of these differ- 
ent control system of arterial pressure is through the feedback gains which 
occur as a result of a given pressure disturbance. Figure 7 illustrates 
a preliminary attempt at such a characterization (Guyton et al (1972b)), 
The solid portions of these curves are based on experimental results, 
while the dashed portions are extrapolations and are somewhat hypothetical. 
Note that the gains fall into three general categories. First, the three 
nervous control mechanisms act very quickly (within seconds on the 


MAXIMUM FEEDBACK GAIN AT OPTIMAL PRESSURE 
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FIGURE 7. 


FEEDBACK GAINS OF SOME PRESSURE 
REGULATORY MECHANISMS. 

FROM GUYTON et al (1972b) 


IS 


essentially logarithmic abscissa). Then several mechanisms (stress- 
relaxation, renln-anglotensln, and fluid shift) begin to act only after a 
few minutes and reach full development within a few hours. Finally, 
the renal-body fluid mechanism and the aldosterone control mechanism 
require hours to days to affect arterial pressure maximally. A remark 
is in order on the nature of the gains quoted In Figure 7. Each of these 
represent a single system result; 1. e. , all systems except the one In 
question are presumably not res|X)nding to the stimulus. Thus, gains 
listed are only Indicative of (not a measure of) expected In vivo responses. 
For example, the Infinite gain characteristic of the renal-body fluid system 
(which implies an integral controller with a physiologically fixed set point) 
only occurs if fluid intake (net drinking) Is absolutely normal and if the 
renal function curve, which relates arterial pressure to urinary output 
rate, remains absolutely normal. If such is not the case, the gain may 
be high, but is net infinite. In spite of the somewhat hjrpothetical nature 
of the gain curves, they are useful as a means of categorizing both the 
relative speed and effectiveness of the pressure regulating mechanisms 
of the body. 

Simulation of more pathological states and some therapeutic 
examples are illustrated in the following. Figures 8-13 illustrate 
typical experiments simulated with this model. These figures show the 
diversity of simulations which this iarge, though simple, model is capable 
of performing. Figure 8 iiiustrates the results obtained by a 2-liter infus- 
ion of isotonic saline into an effectively nephrectomlzed subject (infusion rate 
67 ml/mln. ). Figure 9 pictures the development of hypertension in a 
renal deficient subject during salt loading. Figure 10 shows a simulation 
of Goldblatt hypertension following a renal arterial clamp. Figpire 11 
illustrates the effects of plasma protein deficiency (h 3 rpoproteinemia) on 
circulatory function. Figure 12 shows the simulation of circulatory 
changes following opening and subsequent closing of a massive arterio- 
venous fistula. Figure 13 pictures the development of (bilateral) 



TIME(min.) 


FIGURE 8. ISOTONIC SALINE TRANSFUSION SIMULATION 
FROM WHITE & CROSfON (1974). 
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DAYS 


FIGURE 9. SIMULATION OF HYPERTENSION DEVELOPING 

DURING SALT LOADING OF RENAL DEFICIENT 
SUBJECT. FROM GUYTON et al (1972). 


TOTAL PERIPHERAL RESISTANCE BLOOD VOLUME 
(mm.Hg/liter/min.) (liters) 
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Applied Goldblatt Clamp 



FIGURE 10. SIMULATION OF GOLDBLr.TT HYPERTENSION. 

FROM GUYTON et al (1973b). 



Total Plasma Cardiac 



FIGURE 11. SIMULATION OF THE CIRCULATORY EFFEC 

PROGRESSIVE HYPOPROTEINEMIA. 

FROM GUYTON et al (1972). 
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Open A-V Fistula Close Fistula 
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DAYS 


FIGURE 12. SIMULATION OF CIRCULATORY FUNCTION FOLLOWING 

OPENING AND SUBSEQUENT CLOSING OF A-V 
FISTULA. FROM GUYTON et al (1973). 


Output Pressure Output Pressure Volume 

(ml. /min.) (mm.Hg) (Liters/min.) (mm.Hg) 
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congestive heart failure. Each of these simulations Is realistic In the 
sense that average experimental results compare quite well with the 
representation predicted by the model. 

For the saline Infusion simulation, Figure 8 , the variables shown 
are (top to bottom) plasma volume (VP), interstitial volume (VT8), capil- 
lary pressure (PC), plasma colloid osmotic pressure (PPC), Interstitial 
fluid pressure (PIF), Interstitial colloid osmotic pressure (PTC), capil- 
lary filtration rate (VTC), and lymph flow rate (VTL). The Infusion 
began at time zero and ended at 30 minutes (Input: not shown In Figure 8). 
Note that capillary filtration increases markedly as the transfusion begins 
and that lymph return Increases more slowly resulting In an outward flow 
of fluid from the plasma. At the conclusion of the infusion capillary 
nitration fell rapidly tb equal lymph flow and only 265 ml of fluid was 
transferred to the Interstitial space before equilibrium was obtained. In 
the new steady-state the infused saline was thus distributed in the ratio of 
about 1:3 between the plasma and interstitial space. If the infusion had 
been large enough for the "safety factor for edema" to be exceeded on the 
other hand, almost all of the Infused fluid, beyond a certain point, would 
have entered the Interstitial space. This experiment indicates quite 
clearly the dynamic nature of fluid shift mechanisms and demonstrates 
the suitability of the model for studying such problems. 

Prior to the simulation shown in Figure 9 (about 10 days simula- 
tion), the model subject had its effective renal mass reduced to 30 percent 
of normal (a formulation of analogous pathological state) modifying circula- 
tory variables in the model. A new steady-state developed which was not 
very different from the normal steady-state as long as the diet was kept 
normal. At the break in the curves intake of salt was increased to five 
times normal (Input: not shown in Figure 9), and the simulation pictures 
the classical development of salt-induced, renoprival hypertension. Note 
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In particular the initial rite In cardiac output that leada to an increaalnf 
arterial pressure, In spite of a slight fall Ir. total peripheral resistance, 
and the aubaequent rising peripheral resistance enhancing and then main- 
taining an increased mean pressure in spite of cardiac output returning to 
near normal. Such a simulation illustrates the danger of drawing con- 
clusions about what happens In a situation by examining the steady-state 
only. 

The simulation of Figure 10 illustrates the sequence of events 
following application of a Goldblatt clamp causing sudden renal artery 
constriction in both kidneys (input). Note especially the biphaslc behavior 
of both the cardiac output and peripheral resistance curves. The initial 
large peripheral resistance response is due to the greatly Increased 
angiotensin secretion, while the final state Is primarily due to reduced 
capacitance of the circulation and somewhat increased fluid retention. 

The effects on the circulation of progressive depletion of plasma 
protein (hypoprotelnemia) are shown in Figure 11. At the Initial break 
in the curves the plasma was made to lose protein as shown by the third 
curve from the top (input). 71118 loss decreased the plasma colloid 
osmotic pressure, causing fluid to migrate into the Interstitial space. 

Tlie resultant decrease in plasma volume reduced cardiac output and 
arterial pressure slightly. When interstitial pressure rose to about 
atmospheric pressure (at about the 6-day mark) the Interstitial volume 
began to Increase drastically an;! pitting edema began to develop with 
most of the fluid in the Interstitial spaces being in the free (non-gel) form. 
At about the twelfth day the loss of protein was reduced slightly and this 
had the effect of shifting the capillary flow enough to cause a high urinary 
output and some reduction of the level of edema. 

Figure 12 shows the sequence of events taking place following the 
opening and later closing of a large arterio-venous fistula. At the 
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opening of the fistula (Input), peripheral resistance fell drsmstlcslly to 
about one-half of normal. This was followed by a rise In cardiac output 
(due to enhanced venous return) and a fall In mean pressure. The falling 
arterial pressure elicited an autonomic response which Increased heart 
rate, and the lower renal arterial pressure led to an Increase In renin 
secretion and angiotensin formation. The ultimate effect of the fall In 
pressure then was to decrease urine output and Increase drinking leading 
to an Increase In extracellular (and blood) volume. This Increase 
returned pressure to normal and a steady-state developed with SO percent 
normal peripheral resistance, twice normal cardiac output and normal 
values for most other variables with the exception of fluid volumes. 
Closing of tho fistula (removal of Input) led to essentially opposite effects. 

Figure 13 pictures the effects of repeated reduction In ventricular 
pumping ability of both sides of the heart. Shown are (top to bottom) 
plasma volume (PV), left atrial pressure (LAP), right atrial pressure 
(RAP), cardiac output (CO), free fluid volume In the lungs. (Pul. fluid), 
arterial pressure (A. F.), extracellular fluid volume (E.C.F. V. ), and 
urinary output (U. O^. At the first break In the curves, the pumping 
ability of both ventricles was reduced to 0. 3 nnimal and at each subse- 
quent break in the curves this pumping ability was further reduced by 
about 10 percent of the value just prior to the break (Inputs; not shown In 
Figure 13). After each reduction a marked decrease In cardiac output 
and arterial pressure occurred, but these returned toward normal rapidly. 
Urlnp *7 output decreased and extracellular volume then Increased rapidly 
and the atrial pressures Increased. After a day or so all variables 
tended to return to normal, but with each subsequent attack the transients 
became more severe until finally recovery became insufficient (Note: a 
steady increment of magnitude of responses of these variables, especially 
extracellular fluid volume and a poor response of urinary output In 8-9th 
day episode.) At this point pulmonary edema became severe (due to the 
high left atrial pressure) and death resulted as the oxygen saturation of 
the blood (not shown), decreased below 50 percent. 
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2.2 CURRENT CAPABILITIES 

Examples Lhown here clearly depict the following capabilities 
of existing simulation models; 

1. Tlie ability to simulate analogous pathological states and effects 
of abnormal environmental btrespors by manipulation of system 
variables and changing Inputs In various sequences. 

2. The ability to simulate time courses of responses of controlled 
variables caused by the altered inputs and their relationships. 

3. The ability to simulate physiological responses of treatment such 
as Isotonic saline transfusion as shown in Figure 8. 

‘i. Since the ability of changing Inputs during a course of simulation 
exists, e.g. , opening and closing of A-V fistula in Figure 12, it 
is possible to simulate effectiveness of a treatment, as well as 
to simulate effects of complication superimposed on an existing 
pathological state. 

5. Since sequence and magnitude of input representing a treatment 
can be manipulated, it is pos' ible to compare the effectiveness 
of various treatments/countermeasures for a given pathological 
state. 

Current capability to manipulate input to the whole-body simulation model 
and the supportive data analysis system can have almost unlimited poten- 
tial in its clinical and therapeutic research applications. The use of 
slmulatton models in concert with studies like those mentioned above 
with their ease of manipulation of inputs, the availability of various 
output that is difficult to obtain experimentally, and the computation of 
related system parameters would certainly provide a powerful tool in 
the research of dynamic pathology as well as in the search of therapeu- 
tic means (inputs/countermeasures) and for aid in prognostication of 
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outcome. For example, parameters obtainable from non- invasive 
measurement form only 3-4 percent of total parameters ( 8 para- 
meters out of total 207 physiological important parameters) in Guyton's 
Integrated model. Development of the capability to select and display 
invasive parameters desired, capability of manipulating inputs, and the 
computation of parameters corresponding to both non-lnvasive and 
invasive measurements from subjects when applied to the current 
system would open a way in which models could yield a non-invasive 
.neasurement tool which may be useful for evaluation of patients who 
could not tolerate Invasive procedures. 

In summary, results of this study have exhibited that existing 
simulation models have some capability in clinical applications with 
the few examples presented. However, this work Is in an early stage. 
Broad, sometimes unforeseeable, applications of models to clinical, 
therapeutic, and countermeasure research for NASA's future needs 
represents a new and largely undefined area of research. 
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3.0 


FE/jMUn 9r AFOrTOW T« IMtflATWf mVSkM 

DIAOHOmc AND THIBAPIVTIC RICTAItCH PWOBLIM8 


A preliminary investigation indicates that the feasibility of apply- 
ing current simulation models in diagnostic research requires further 
basic study. Other types of models, such as stochastic models, 
probably have greater applicability to clinical diagnosis. However, 
since the contract work to date has considered only simulation models 
of major body subsystems, this report considers only applications of 
this type of model. A few areas of possible application are noteworthy 
because of their theoretical applicability even if their practicality 
remains in question. These areas include utilization of sensitivity 
coefficient analysis methods, particularly stability analysis, and 
inyerse sensitiyity analysis using current simulation models. Appli- 
cations in therapeutic research appear to proyide the most fruitful 
area for future studies and some of these ideas are discussed in this 
section. 


3. 1 SIMULATION MODEL APPUCATIONS IN CLINICAL DIAGNOSIS 

Analysis of the stability of biological systems becomes a useful 
aid to clinical diagnosis if stability can be measured and recognized 
quantitatiyely for major body subsystems. Examples of such needs 
are assessment of crew health status and recognition of a change in 
pathological state (e. g. , recorded as the daily progress report of a 
patient in clinical medicine). 

Theoretically, the stability of a dynamic system has an inyerse 
relationship with sensitivity in a negative feedback system. In general, 
sensitivity to disturbing factors can be reduced by an Increase in feed- 
back gain. On the other band, the onset of instability occurs as a 
consequence of this gain increase. Thus, systems with high gain may 
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have low aeneitivlty to external perturbations. Most biological 
systems are normally stable and they exhibit low sensitivity. Whether 
or not they are working somewhere near the stability limit via high 
gain factors is not known, but should be utudled on a case-by-oase 
basis (Jones, 1973). Inherent Instability is dependent upon the 
properties of the system and Is mx normally a function of the specific 
disturbance. If the system is Inherently stable, all transients will 
ultimately disappear regardless of the disturbance causing them 
although biological systems have their own threshold to the distur- 
bance. On the other hand, any disturbance to an unstable system 
will initiate oscillations that Increase in amplitude with time. Insta- 
bility can arise from either Inherent features of the real system or 
from structural features of the mathematical model. Techniques of 
recognizing such instability of selected systems appears feasible util- 
izing sensitivity analysis methods which can reveal both types, those 
inherent to the real body system and those from the simulation model 
Itself although It is not always possible to distinguish between the two. 

In clinical diagnosis, whether the biological system is stable or 
unstable may not be important, but the degree of stability and its 
change may be the concern. An analysis of stability, thus, not only 
becomes an important measure of the competence of a model, but £.iso 
becomes a useful aid in systematic diagnostic research when there is 
confidence in the mathematical simulation. Little work has been done 
on stability analysis of complex physiological systems ^ such a prac- 
tical nature. Formal techniques for investigating stability in linear 
systems and for simple nonlinear systems have been reported (Jones, 
1973; Rosen. 1970; Tomovlc, 1963; and Tomovlc, et al , 1972), but for 
the most part, studying large scale nonlinear models is a trial and error 
experience. Applications of the simulation model in this area, there- 
fore, requires further study. 
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Inverse sensitivity analysis is another technique which appears 
to have a potential application in diagnostic research. If the problem 
of sensitivity analysis is expressed as determining the behavior of a 
model, given all the parameter variations, then the Inverse problem 
would be to determine (or identify) the parameter variations capable 
of producing a given behavior of the real system. The problem of 
direct sensitivity can be solved satisfactorily if the sensitivity coeffic- 
ients of the dynamic system are known. The method for finding 
inverse sensitivity is also based on a knowledge of the sensitivity 
coefficients. Unlike direct sensitivity analysis, inverse sensitivity 
analysis requires data measurements from the real s}rstem. 

In essence, solving the inverse sensitivity problem is not unlike 
parameter estimation analysis. Analytic methods have been worked 
out for the linear system (Tomovic, 1963) while empirical methods 
must be used for nonlinear systems. With Inverse senslt! / tty, one 
starts from the known perturbation of the transient state and seeks the 
set of parameter perturbations capable of causing them. This inverse 
problem may not have a unique solution. Uniqueness is dependent 
primarily on the number of parameters vs. the number of variables 
measured and the number of points in time at which measurements 
are made. Nevertheless, it would be valuable to know the various 
solutions possible since this would narrow down the range of search 
for the possible causes of observed pathological or undesirable states. 

If several different parameter perturbations could produce similar 
model results, it may be possible to accept the most reasonable baseo 
on physiological plausibility or alternatively this Information could 
provide the basis for further experiments and diagnostic testings. 

It should be noted that the clinical diagnostic procedures do not 
necessarily follow the procedure of a direct determination of Inputs 
(causes) from the observed responses in the strict sense of a control 
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system. The clinlcsl definition of "diagnosis" is a<fo rather vague. 
Physicians use a treatment as a means of diagnosis in some cases. 
Inverse sensitivity analysis could, therefore, provide a more syste- 
matic approach to the diagnostic processes. 

It is worth mentioning here that numerous studies have been 
reported on computer-aided medical diagnosis which are based on 
stochastic approaches. This area of investigation was not included 
as a part of this contract. A simple stochastic compterizod diagnosis 
assistance program based on a diseases-symptom matrix to investigate 
and demonstrate the feasibility of this approach has been developed as 
a part of an independent research and development project. Th s pro- 
gram is currently capable of making differential diagnosis of 336 
diseases from among 82 signs and symptoms. 

3. 2 THERAPEUTIC/COUNTERMEASURE RESEARCH 

Applications of current simulation models in therapeutic/counter- 
measure research is feasible and may have applicability in broad areas 
of interest. An example of such interest was indicated in the recent 
Skylab Life Sciences Symposium. The need for developing counter- 
measures which were indicated by the Symposium include: 

o Countermeasures for habitation space sickness in the 
first 5 days of space flight. 

o Countermeasures for fluid redistribution or orthostatic 

intolerance (such as hypotensive /hypertensive suit) for 

Shuttle passenger population groups. 

o Countermeasures for cardiovascular deconditioning. 

o Countermeasures for the maintenance of muscular mass 
and functional integrity against their loss in space flight. 
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o Countermeaaures for inflight bone calcium Iota, 
o Countermeasures for the erythropoietic inhibition/ari ant 
in space. 

o Countermeasures for body weight loss in space. 


Many nreas of therapeutic research could be Investigated with little or 
no modification* to the current models. This section Illustrates the 
feasibility of this application with a few examples. 

D. C. Riggs (1970) describes therapeutic ta~‘ge^ points in his 
explanation of vicious circles using a case of paroxysmal nocturnal 
dyspnea di£ to acute pulmonary edema with a symbol-and-arrow diagram 
as shown in Figure 14. When a vicious circle underlies a pathological 
process, measurements of the open-loop gain would warn of impending 
disaster and would serve as a quantitative index of the effectiveness of 
treatment. Unfortunately, the variables in the positive feedback loop 
are often difficult to identify and even more difficult to measure accura- 
tely enough so that feedback equations describing the interrelations 
between two members of the loop could be written. 

With left heart failure, acute episodes of pulmonary edema charac- 
teristically occur at night after the patient has gone to bed. Unless 


prompt treatment is given, a severe attack may be fatal presumably 
because the overall gain of the vicious circles equals or exceeds unity. 

In Figure 14, the steps directly involved in the pathogenesis of 
pulmonary edema and the consequent anoxia have been identified by 


heavy arrows. In left heart failure, resistance of the 

left ventricle to the inflow of blood is increased, so that ^Latr* 
pressure in the left atrium is abndrmally elevated, necessitating an 


elevation of P, and hence of P,.,,, the net transcapillary 

lung capill flit 

filtration pressure in the lungs, which causes movement of fluid from 


blood plasma into the alveoli. During the day, while the patient i<* 
"up and around", the force of gravity tends to limit venous return. 
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R • Resistance of the left ventricle to the inflow of blood 

L venf 

D 

periph: Peripheral resistance 

p 

L atr: Pressure in the left atrium 

p 

Lung capill; Pressure in the lung capillaries 

'^R vent: Output of the right heait (venous return) 

P^.|^: Net transcapillary filtration pressure in the lung 

p 

air flow: Resistance of air flow in the lungs 
Vaiv : Alveolar ventilation 

Pq 2 art* pressure of O 2 in arterial blood 

Yfilt * transudation into the alveoli 

^002, art- pressure of CO 2 in arterial blood 

W 

I’csp; Rate of respiratory muscle work 
%2 cons* ^2 

9qo 2 prod* ^^2 


FIGURE 14. ACUTE PULMONARY EDEMA (RIGGS, 1970) 
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(V{^veQt* output of the right heart) by pooling fluid In dependent 

parta of the body, and by vaaoconatrlctlon (mediated via the baro- 

receptora) which Increaaea peripheral realatance, R , . . (Veno- 

peripn. 

motion, however, tends to maintain venous return by decreasing the 
volume of blood pooled In the veins). When the patient goes to bed, 
the angle which the long axis of the body makes with the horizontal la 
decreased towards zero; gravitational forces no longer hinder the 
return of blood to the heart, and cardiac output Increases, particularly 
If there Is any dependent edema which can be reabsorbed Into the 


Pi .11 ®nd 

lung caplll flit 


As Increases, the rate of fluid 


transudation Into the alveoli becomes more rapid and the volume of 
edema fluid In the lungs Increases. Fluid accumulation increases 
the resistance to air flow into the alveoli (thus reducing the 

alveolar ventilation), and decreases the diffusing capacity of the lungs. 
Both effects reduce p^^ partial pressure of oxygen In arterial 

blood, and cause anoxia. Both effects also Increase p 


C02,art 


, the 


partial pressure of carbon dioxide In arterial blood. 

Key portions of positive feedback loops are emohasized by thick, 
but unfilled arrows. All of these positive loops depend upon the ability 
of hypercapnia and hypoxia to Increase W , the rate at which the 
muscles of respiration do work. The positive feedbaci' loops which 
Involve Increased O consumption and CO production by the respiratory 
muscles were described In his model on the chemical control of breath- 
ing. In pulmonary edema, two other positive feedback paths establish 
additional vicious circles. First, the labored, gasping respiration 
(dyspnea) decreases the mean Intrathorarlc pressure. This decrease 
not only tends to promote return of blood to the heart, but also decreases 
the pressure within the alveoli, thereby Increasing Second, as Is 

true of any kind of exercise, the Increased work of the respiratory 
muscles will cause an Increase In cardiac output via mechanisms which 
are still not entirely clear. 
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The preponderance of poaltlve feedback in pulmonary edema caueea 
a paradox. The rapidly worsening tissue anoxia would seem to demand 
measures to increase delivery of oxygen from the air to the lungs via 
Increased lung ventilation, and from the lungs to the tissues via Increased 
cardiac output; yet Increased cardiac output and increased ventilatory 
efforts are two of the main factors accounting for the initiation and 
perpetuation of the pulmonary edema. In fact, as clinicians have long 
recognized, in pulmonary edema It is essrntial to interrupt the greater- 
than-unlty-gain positive feedback loops even if we do so by depressing 
respiration and decreasing cardiac output. 

The following list shows what therapeutic measures are available 
to reverse the rapidly progressive deterioration characteristic of 
untreated pulmonary edema. Each circled number in Figure 14 shows 
the point 

o To 

1 . 


2 . 

o To 

3. 

o To 

4. 

5. 


of attack of the correspondingly numbered therapeutic measure, 
interrupt the vicious circles 

Morphine sulfate Intravenously. Lowers by decreas- 

ing the sensitivity of the respiratory center to stimuli caused 
by hypercapnia and hypoxia. 

Assisted respiration. (Gas mixture is supplied at greater- 
than-atmospheric pressure during inspiration.) In part, 
substitutes exogenous for endogenous respiration work. Also 
taises mean intrathoracic pressure. 

combat anoxia directly. 

Increasing the partial pressure of oj^gen in the inspired air. 

decrease venous return and cardiac output 
Bleeding (venesection). 

Pooling blood in limbs by applying tourniquets to obstruct 
venous outflow ("bloodless venesection"). 
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6. Elevation of head of bed, or altting poeture, to promote 
gravitational pooling. Thla may alao favor redlatrlbutlon 
of fluid In the lungs ao that at least the upper portions can 
be aerated, thus Increasing diffusing capacity and decreaa- 
Ing reaistance to air flow. 

o To combat the heart failure. 

7. Amlnophylllne slowly Intravenously, for Its positive Inotropic 

effect on the heart, and to relax bronchlolar smooth muscle 

thus reducing R , » 

® air flow 

8. Digitalization with a rapidly-acting cardiac glycoalde, to 
relieve the heart failure. 

o Longer term preventive measures 

9. Maintain digitalization with a suitable drug. 

10. Administer diuretics If needed to prevent accumulation of 
dependent edema. 

The above example suggests a method of applying a simulation 
model In therapeutic/countermeasure research which Is: (1) the 

requirement to interpret effects of physical and chemical (Including 
drugs) treatments to their biological effects which can be expressed 
by variables In the model, (2) their target points have to be identified, 
and (3) a simulation model for the above requires a program of the 
vicious circle with a regulation threshold model which controls Inclusion 
or exclusion of the vicious circle. 

Recently, Interest has developed concerning whether or not the 
existing hypertensive anti-g suit successfully used for Skylab missions 
is adequate for the Shuttle passenger population. Clinical applications 
of hypertensive and hypotensive devices applied to the lower body have 
been reported by many researchers (e.g. , by Wothuis, Bergman, and 
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Nicogossian, 1974). In addition to static hypertensive or hypotensive 
lower body devices, an external counterpulsation device (ECP) has 
been applied to combat against cardiogenic shock of myocardial infarc- 
tion, victims. (Birtwell, et al, 1969; Soroff, et al, 1969, 1971, and 
1974). Counterpulsation raises the diastolic pressure end lowers 
the systolic pressure to assist circulation. It accomplishes this by 
synchronizing its pumping action with the heart's own action triggered 
by the R wave of electrocardiogram. The pumping pressure is 
exerted in the opposite way creating blood pressure of 80 over 120 mm Hg. 
The hemodynamic response should be noticeable within 45 minutes after 
the use of positive-negative pressure assist. By raising the diastolic 
pressure, myocardial and coronary perfusion is increased. The work 
of the heart is reduced by taking the pressure off the left ventricle 
thereby increasing its stroke volume. Improved perfusion oxygenation 
of the injured myoc irdium and providing a rest to the impaired left 
ventricle are the principal therapeutic aims of this approach. Soroff, 
et al (1969) also claims the evidence of the development of more and 
larger collateral myocardial vascularization. The authors report 
that there is no evidence of impairment of peripheral circulation or 
aggravation hemolysis and no reason to anticipate thromboembolism. 

The only contraindication for its use is in patients with aortic insuffi- 
ciency. 

The cardiovascular model for LBNP simulation is capable of 
simulting the effects of such devices with only minor changes. The 
evaluation of existing and proposed devices for Shuttle using the simu- 
lation model will not only reveal effects of such devices, but should 
also indicate the most effective approach before actual construction is 
required. This evaluation should consider inflight prophylactic 
devices as well as postflight therapeutic devices. 
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4.0 and RECOMMENDATMMl 

ITie feasibility of applying subsystem and whole-body simula- 
tion models to diagnostic and therapeutic considerations has been 
studied, and examples of model applications are given In Section 2.0. 

It Is concluded that simulation models can bo applied to clinical and 
aerospace medicine areas, although this work Is In a very early and 
preliminary stage. The results here indicate that there are many 
areas to be explored and models to be developed to produce practical 
applications to diagnostic and therapeutic problems in clinical research 
as well as aerospace medical operations. 

Baaed on the preliminary otudy in this area, the following study 
tasks are recommended for a tnore in-depth and advanced research 
plan; 

(a) Select applicable pathological states, or diseases, and 
injuries of Interest to N/SA aerospace operations. 

(b) Survey and formulate known physiopathological mechanisms 
and Interpretation of treatments into variables and target points for 
development of pathological simulation models including vicious circle 
models. 

(c) Perform a system requirements analysis for the applica- 
tion of diagnostic and pathological simulation models to an integrated 
crew health monitoring system for the Shuttle era. 

(d) Investigate the possible use of stochastic models in com- 
bination with the whole-body algorithm to provide a total system diag- 
nostic and therapeutic model. 

(e) Develop a model to simulate effects of lower body counter- 
pressure devices. 
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(0 Perform ■ImulatloM of constant LBNP effects on body 
fluid redistribution after orbital insertion. 

(g) Survey past studies and investigate simulation of constant 
pressure hypertensive g-suit used during reentry and post-landing 
phases. 
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